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ABSTRACT 

We present new observations of the plerionic supernova remnant 3C58 with the VLA at 74 
and 327 MHz. In addition, we re-reduced earlier observations at 1.4 and 4.9 GHz taken in 1973 
and 1984. Comparing these various images, we find that: 1. the remnant has a flat and relatively 
uniform spectral index distribution, 2. any expansion of the remnant with time is significantly less 
than that expected for uniform, undecelerated expansion since the generally accepted explosion 
date in 1181 A.D., and 3. there is no evidence for a non-thermal synchrotron emission shell 
generated by a supernova shock wave, with any such emission having a surface brightness of 
< 1 x 1CT 21 W mr 2 Hz" 1 sr" 1 at 327 MHz. 

Subject headings: ISM: individual (3C58) — supernova remnants 



1. INTRODUCTION 

Plerionic, or filled center supernova remnants 
(Weiler & Panagia 1978; Weiler & Sramek 1988) 
have become established as a class. The supernova 
remnant 3C58 (G130. 7+3.1) was an early mem- 
ber of the class, and it has often been described as 
morphologically similar to the Crab Nebula, the 
Crab being the prototype for this class. 3C58 is 
at a similar distance as the Crab, 3.2 kpc (Roberts 
et al. 1993), has a similar size, 6' x 9', and shows 
a similar, filamentary, filled center morphology 
in the radio (Weiler & Seielstad 1971; Wilson & 
Weiler 1976; Green 1986; Reynolds & Aller 1988). 
Like the Crab, 3C58 has a relatively flat radio 
spectrum with a = —0.10 ± 0.02 (S oc v +a ) from 
40 MHz to 15 GHz (Wilson & Weiler, 1976; Green 
1986). 3C58 has been moderately well identified 
with the historical supernova event of 1181 A.D. 
(Clark & Stephenson 1977), which would give it an 
age similar to that of the Crab, which is associated 



with SN1054. 

Plerions are thought to be powered by pulsars, 
which build the non-thermal nebula through the 
injection of relativistic particles and amplification 
of magnetic field, although the exact injection and 
acceleration processes are not yet well understood. 
In the optical, radial velocities < 900 km s -1 have 
been measured for the filaments in 3C58 (Fesen 
1983), which is consistent with the association 
with SN1181 only if the remnant has been sub- 
stantially decelerated or the filaments are due to 
the excitation of in-situ matter overrun by the su- 
pernova shock. However, the linear size of < 9 pc 
(for a distance of 3.2 kpc) along with the age of 
only 820 years imply rather low expansion veloci- 
ties, which suggests a low energy supernova event. 
In addition, the radio flux density has been shown 
to be increasing (Aller & Reynolds 1985; Aller, 
Aller, k Reynolds 1986; Green 1987). It seems 
(Salvati et al. 1998; Woltjer et al. 1997) that a 
non-standard evolution of the pulsar's output is 
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required to reconcile this and the relatively low 
frequency of the observed spectral break in the 
synchrotron emission (Green & Schcuer 1992) with 
an age of 820 years. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

We observed 3C58 using the A and B configu- 
rations of NRAO's Very Large Array (VLA) at 
74 MHz using 63 spectral channels for a total 
bandwidth of 1.5 MHz, and at 327 MHz using 31 
spectral channels for a total bandwidth of 3 MHz. 
Spectral line mode was used to be able to excise 
radio frequency interference (RFI). 

The data were calibrated and reduced using 
NRAO's AIPS software package. The instrumen- 
tal bandpass was determined using observations 
of Cygnus A. The amplitude gains were set using 
Cygnus A at 74, and 3C48 at 327 MHz, and tied 
to the flux scale of Baars et al. (1977). Reduction 
of 74 MHz data requires somewhat different tech- 
niques (see Kassim, Perley, & Lazio 2000). RFI is 
a significant problem and the data were first auto- 
matically edited for deviant samples in frequency 
space by the AIPS task SPFLG. We then averaged 
the data into a manageable number of channels 
and further edited each channel manually. 

Sidelobes from sources well outside the pri- 
mary beam are a significant problem at 74 MHz 
and Cassiopeia A, which was 20° away, pro- 
duced strong sidelobes in our field-of-view. In 
order to correctly remove these sidelobes, three- 
dimensional imaging is required (see Perley 1999) 
and the AIPS task IMAGR was used to simulta- 
neously image and deconvolve numerous flat but 
non-coplanar regions. These were then reassem- 
bled into a single image by the task FLATN. 

3. RESULTS 
3.1. Images 

The full resolution image at 74 MHz is shown in 
Figure 1 . This image was made from the combined 
A and B array data sets with uniform weighting, 
and had a restoring beam size of 26" (FWHM). 
The peak brightness for 3C58 was 0.55 Jy bm _1 , 
the background rms was 20 mJy bm _1 , and the 
total flux density for 3C58 was 33.6 Jy. The im- 
age at 327 MHz is shown in Figure 2. It was 



made from the combined A and B array data sets, 
and had a restoring beam size of 8.2" (FWHM), a 
peak brightness for 3C58 of 63 mJy bm _1 , a back- 
ground rms of 0.57 mJy bin -1 , and a total flux 
density for 3C58 of 33.9 Jy. It should be noted that 
no compact emission source is detected near the 
center of the nebula. In addition, in Figure 3, we 
show a low resolution image of the whole primary 
beam area at 327 MHz, convolved to a resolution 
of 3.3'. No shell emission is visible. The average 
background rms between the edge of the field and 
the center is 37 mJy bm _1 , or 3.4 x 10~ 22 W m~ 2 
Hz -1 sr _1 . We can accordingly place a 3cr limit on 
any shell emission of 1 x 10~ 21 W m~ 2 Hz -1 sr _1 
at 327 MHz. While the field of view at 74 MHz is 
even larger, the sensitivity is considerably poorer, 
so the higher frequency gives our most sensitive 
limit on extended emission. We note however, 
that the 74 MHz observations would be sensitive 
to very large shells, i.e. with diameters > 1.5° and 
we see no such shell emission. 

3.2. Spectral Index Distribution 

We made an image of the spectral index be- 
tween 74 MHz and 327 MHz 1 , ag;j$, as follows. 
The 327 MHz image was convolved to the larger 
resolution at 74 MHz and then the two images 
were aligned by the peaks of the point sources 
near 3C58 (at 74 MHz, ionospheric refraction can 
cause substantial displacement of the entire im- 
age). This image of o^of, which has a resolu- 
tion of 26", is shown in Figure 4 (regions where 
the uncertainty in o^of was larger than 0.4 have 
been blanked for clarity). The spectral index has 
a mean value of —0.04 over the nebula, with an 
rms of 0.18. Only random variations are visible in 
Figure 4. In particular, no steepening towards the 
edge of the nebula is observed. 

We can calculate the integrated value of a^of 
from the total flux densities in the naturally 
weighted images, and we obtained o^gy (integrated) 
= —0.07 ± 0.05, where the uncertainty comes 
chiefly from the uncertainty in the VLA ampli- 
tude calibration. This value is consistent with the 
"o 5 04 = -°- 10 ± °- 02 S iven b y Green (1986). 

We also re-edited and re-imaged the data of 
Reynolds & Aller (1988) and made a spectral in- 

1 We will refer to spectral indices with sub- and superscripts 
indicating the relevant frequencies in GHz. 
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dex image between 327 MHz and 4.9 GHz, 0^33, 
which is shown in Figure 5. Again no particular 
features are visible in the body of the nebula, and 
no steepening towards the edge is observed. The 
mean and rms scatter of the spectral index over 
the central region of 3C58 were —0.051 and 0.014, 
respectively. There is a very slight suggestion of 
a steepening of the spectral index near the fila- 
ments by < 0.02, which is likely not real, since 
the systematic sources of error, in particular im- 
age fidelity limits (see e.g., Perley 1999), will in- 
troduce errors in the spectral index at this level 
which are likely correlated with the structure in 
the image itself. The integrated value of 0:^33 was 
—0.06 ±0.03, consistent with the above value from 
Green (1986). 

4. EXPANSION OF 3C58 

In order to determine the rate of expansion of 
the synchrotron nebula, we compared our results 
to those from the aforementioned observations in 
1984 by Reynolds & Aller (1988), and from obser- 
vations in 1972 by Wilson & Weilcr (1976). We 
re-edited and re-imaged both these sets of obser- 
vations using MEM deconvolution. In Table 1 we 
give details of the images used for determining the 
expansion. All of the images at 1.4 and 5 GHz have 
been corrected for the attenuation of the respec- 
tive primary beam patterns (at lower frequencies 
the correction is negligible). Our goal was to de- 
termine the expansion speed of the radio nebula. 
Since there are few well-defined, compact features, 
we choose to measure the expansion not by deter- 
mining the proper motion of individual features, 
but by determining an overall statistical scaling 
between our images. We accomplished this by us- 
ing the MIRIAD task IMDIFF which calculates 
how to make one image most closely resemble an- 
other by calculating unbiased estimators for the 
scaling in size, e, the scaling and the offset in 
flux density, A and b respectively, and the offsets 
in RA and Dec, x and y respectively, needed to 
make the second image most closely resemble the 
first. We are chiefly interested in the expansion 
factor, e, but because of uncertainties in flux cal- 
ibration, absolute position, and image zero-point 
offsets caused by missing short spacings, we deter- 
mined all parameters. This method was developed 
by Tan & Gull (1985) and more details of its use 
in a similar situation are given by Bietenholz et 



al. (1991). 

Since we have shown above that the radio spec- 
tral index is uniform to within the uncertainties, 
images at any two frequencies, v\ and v 2l should 
closely resemble one another, except for a flux scal- 
ing factor which is given by (^1/^2)", where a is 
the spectral index. Thus we can usefully compare 
images taken at different observing frequencies to 
test for expansion. 

We ran IMDIFF on various combinations of the 
above images, always with both images convolved 
to the same effective resolution, to get the expan- 
sion factors shown in Table 2. In order to estimate 
the uncertainty in our calculation of e, we also ran 
IMDIFF on two pairs of (almost) contemporane- 
ous images, for which e should be 0.0. They are 
the expansion between 1.4 and 4.9 GHz at essen- 
tially the same epoch, and the expansion between 
a CLEAN and a MEM image of the 1973 4.9 GHz 
data, and the results are shown in the last two 
lines in Table 2. In both these test cases we ob- 
tain expansions of < 0.2%, which suggests that 
the uncertainties in our expansion factors are of 
this order. 

In order to determine the expansion speed, we 
fit a straight line through the first seven expansion 
factors in Table 2 by weighted least squares. The 
expansion for zero time difference must be 0.0, so 
we constrain the fit line to pass through the ori- 
gin. Since we do not have formal uncertainties 
in the expansion factor, we use the weights listed 
in Table 2, which we derived as follows: it seems 
reasonable to suppose that the uncertainty in the 
expansion factor is proportional to the rms of the 
residuals to the IMDIFF fit, am, divided by the 
average brightness in the images - more flux in 
the image will allow a relatively better fit for the 
same rms. Since the total flux density is roughly 
the same in all the images, the average brightness 
will be inversely proportional to the number of 
beam areas over the source. It also seems reason- 
able that the uncertainty is inversely proportional 
to the square root of the number of beam areas in 
the fitting region, Acting region- We thus adopt as 
the weight 



^fitting region 

For the 5 GHz, 1973 WSRT data we use only ex- 
pansions derived from a version of the image high- 
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pass filtered with a Gaussian of 1' FWHM, be- 
cause this image had evident short-spacing prob- 
lems (Wilson & Weiler 1976). The points, as well 
as the weighted least squares fit, are plotted in 
Figure 6. 

The straight line fit to the data in Table 2 by 
weighted least squares is plotted as the solid line in 
Fig. 6. The slope, or average expansion speed, we 
obtain from the fit is of 0.020±0.008% yr' 1 , which 
corresponds to - 910 ± 360 and 550 ± 220 km s" 1 
along the major and minor axes of the nebula, 
respectively. We can estimate the uncertainties 
from the scatter of the data points about the fit 
line, and we find that a W = 1 corresponds to 
an uncertainty of 0.34% in the expansion factor. 
Our two check points, which measure the expan- 
sion between almost simultaneous images (the last 
two lines of Table 2), are indeed consistent with 
an expansion of 0.0, as they should be. We also 
note that an unweighted fit gives almost the same 
expansion speed and uncertainty. 

Our measured expansion speed (solid line in 
Fig. 6) is much less than would be expected for 
undecelerated expansion of 3C58 since 1181 A.D. 
The undecelerated "historical" expansion speed 
would be 0.124% yr _1 and is plotted as the dashed 
line in Fig. 6. In addition, in Figure 7, we show 
profiles drawn through the long axis of 3C58 at 
1.4 GHz in 1984 and at 0.3 GHz in 1998, and for 
comparison a hypothetical profile from undeceler- 
ated expansion of the 1984 profile to 1998. This 
last can be seen to be notably larger than the true 
1998 profile, also indicating a lower than "histori- 
cal" expansion speed. 

The nebular structure is dominated by the rel- 
atively sharp outer edges rather than by any well 
defined filaments or clumps, so our expansion re- 
sults derive in considerable part from the overall 
extent of the nebula. Unfortunately, this makes 
our expansion results sensitive to short-spacing 
problems and also to the accuracy of the beam cor- 
rections, especially at 5 GHz, where the VLA and 
WSRT primary beams are 825 and 650" FWHM, 
respectively. Because of this we cannot exclude 
the presence of systematic biases in our results. 
However, we feel that the aggregate expansion 
speed, which is derived from a variety of combi- 
nations of different epochs, resolutions, telescopes, 
and frequencies is unlikely to be greatly biased. 



5. DISCUSSION 

We measure an expansion rate of the radio neb- 
ula for 3C58 of 0.020 ± 0.008% yr" 1 (correspond- 
ing to - 910 ± 360 and 550 ± 220 km s^ 1 along 
the major and minor axes of the nebula, respec- 
tively). This is much smaller than would be ex- 
pected if 3C58 were indeed the undecelerated rem- 
nant of SN1181, which would require an expansion 
of 0.124% yr- 1 (~ 5500 and - 3300 km s" 1 along 
the major and minor axes respectively). Con- 
versely, our measured expansion rate suggests an 
age for 3C58 of 5000 yrs with a 3er lower limit of 
2250 yrs, which is incompatible with the generally 
accepted association with SN1181. Although, we 
cannot entirely rule out systematic sources of error 
in our determination of the expansion, we think it 
unlikely that they are large enough to make our re- 
sults compatible with an expansion age of ~ 820 
years. 

This low expansion speed is corroborated by the 
proper motions of optical filaments: Fesen, Kirsh- 
ner, & Becker (1988) and van den Bergh (1990) 
looked for such proper motions and failed to find 
any larger than 0.07" yr" 1 . Fesen (1983) mea- 
sured optical radial velocities for several filaments 
in 3C58. The largest radial velocities he measured 
were ~ 900 km s _1 , for filaments situated near the 
projected center of the nebula, where the expan- 
sion is along the line of sight, and thus where the 
expansion speed is expected to be approximately 
equal to the radial velocity. This is also substan- 
tially lower than the > 3000 km s _1 implied by 
undecelerated expansion and a distance of 3.2 kpc. 
(We note that Fesen, who used a somewhat smaller 
distance to 3C58 of 2.6 kpc, reconciled these ra- 
dial velocities with an age of 820 years by suggest- 
ing that substantial deceleration had taken place, 
and that the average speed since 1181 A.D. was 
~ 3000 km s _1 ). In order to determine an ex- 
pansion age from this radial velocity, we need to 
know the radius at which the measured filament 
lies from the center of the nebula. If we take this 
radius to be the mean projected radius of ~ 4 pc, 
these radial velocities also imply a large age of 
~ 4400 years, but since the optical emission is 
faint, those radial velocities may have been mea- 
sured for filaments nearer the center of the neb- 
ula, and thus a smaller radius and consequently 
a smaller age is possible. However, we consider it 
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unlikely that the largest radial velocities were mea- 
sured for filaments only ~ 1 pc from the center of 
the nebula. If 3C58 is associated with SN1181, 
all this evidence suggests that it must have been 
greatly decelerated. 

The simplest solution to these inconsistencies 
is to suggest that 3C58 is not the remnant of 
SN1181. This also alleviates the difficulties in ac- 
counting for the very sharp and rather low fre- 
quency spectral break observed in 3C58 (Woltjcr 
ct al. 1997; Green & Scheuer 1992). Note that a 
lower value for the somewhat uncertain distance 
estimate for 3C58 would reduce the age deter- 
mined from the radial velocities, but not that de- 
termined from the proper motions, so that is not 
a solution. The positional coincidence between 
SN1181 and 3C58, however, was confirmed by a re- 
cent re-examination (Stephenson & Green 1999), 
and so if 3C58 is not to be associated with SN1181, 
it must be asked where the remnant of SN1181 is, 
that there are no other obvious candidates, and 
that it seems unlikely that the remnant would al- 
ready be too faint to observe at an age of only 820 
years. 

Van den Bergh (1990) suggests that the presently 
visible optical filaments predate the supernova 
event, and are quasi-stationary but energized by 
the passing supernova shock wave. It is not clear 
that the low expansion rate of the synchrotron 
bubble can be accounted for in this way, since it 
seems unlikely that these pre-existing filaments, 
which have a low filling factor, would be able 
to confine the synchrotron bubble. While some 
coupling between the synchrotron bubble and the 
thermal filaments would be expected, since the 
magnetic field in the bubble would wrap around 
the thermal filaments, one would expect the syn- 
chrotron bubble to rapidly expand beyond the fila- 
ments along with the rest of the supernova ejecta. 
The outside edge of the visible synchrotron bub- 
ble should be expanding at least at the speed of 
the local ejecta, which would be ~ 4000 km s _1 
for undecelerated expansion. Should pre-existing 
filaments nonetheless be somehow confining the 
synchrotron bubble, then the filaments should be 
accelerated. The current measurements, unfor- 
tunately, do not have the precision necessary to 
determine if this is the case. 

A possible explanation for slow-moving fila- 
ments is the following. The magnetic field will 



wrap around the pre-existing filaments, and cause 
an intensification of the synchrotron emissivity in 
their vicinity, which is why these filaments are vis- 
ible in radio synchrotron emission. If the structure 
of the radio nebula is dominated by these irregu- 
larities, then our determination of the expansion 
would determine the expansion precisely of the fil- 
aments and not of the synchrotron bubble as a 
whole. We note that expansion speeds we measure 
at lower effective resolutions tend to be slightly 
smaller than those at larger resolutions, as would 
be expected in this case, but even the measured 
expansions at lower resolutions are far below the 
expected values for undecelerated expansion. Fur- 
thermore, profiles through the images do not sup- 
port the conjecture of quasi-stationary filaments 
in an expanding bubble: the profiles displayed in 
Figure 7 show that even the expansion of the edge 
of the nebula appears to be considerably less than 
is expected from undecelerated expansion. 

Let us turn to some general physical consider- 
ations about the synchrotron bubble. The pres- 
sure in the synchrotron bubble is relatively high: 
standard arguments give a minimum pressure of 
> 10~ 10 dyn cm~ 2 . This is the minimum pressure 
of the synchrotron fluid assuming a filling factor of 
1.0, and a line of sight depth equal to the long axis 
of the projected nebula. Both a lower filling factor 
and the admixture of any thermal material would 
raise the pressure, and we think it unlikely that 
the line of sight dimension is much larger than the 
long axis of the projected nebula. 

The slow expansion of the bubble therefore ne- 
cessitates some confining pressure, since the sound 
and Alfven velocities in the bubble are also high. 
This leads to the question of what surrounds the 
presently visible nebula. The three possibilities 
are a) the ISM b) the wind bubble of the pre- 
supernova progenitor of SN1181, or c) a freely 
expanding stellar envelope thrown off in the su- 
pernova explosion. This last might be expected, 
since the total energy in the visible remnant is 
far below the expected energy of 10 51 ergs re- 
leased in a supernova for any reasonable mass in 
the filaments. This last situation does seem to 
obtain in the Crab Nebula, where strong circum- 
stantial evidence points to the existence of such 
a freely expanding stellar envelope (Hester et al. 
1996; Sankrit & Hester 1997), even though the ra- 
dio emission from this shell has so far eluded de- 
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tection (Bictenholz et al. 1997; Frail et al. 1995) 
and the direct optical detection is controversial. 

The static pressure of any of these three ex- 
terior media, p cxt is expected to be considerably 
smaller than the minimum pressure in the syn- 
chrotron bubble: 

p cxt = 1.4xlCT 12 dyn cm" 2 [n/lcm^ 3 ] [T/10 4 K] 

(2) 

where n is the number density and T is the tem- 
perature. This suggests that the synchrotron bub- 
ble is confined by ram pressure as it expands into 
the surrounding material. We are then brought 
to a crucial difference between 3C58 and the Crab 
Nebula: in the case of the Crab, the synchrotron 
bubble has been accelerated since the supernova 
event (Bietenholz et al. 1991), and is therefore ex- 
panding into the stellar envelope, and confined 
by the ram pressure of this expansion. In the 
case of 3C58, by contrast, the synchrotron bubble 
is expanding considerably more slowly than any 
freely expanding material thrown off in SN1181, 
and cannot therefore be confined by ram pressure 
in any freely expanding ejecta. 

This implies for 3C58 that there is no freely 
expanding envelope, making it a very low energy 
supernova event, or that the envelope has been 
greatly decelerated, which seems unlikely in light 
of the non-detection of any shell emission exterior 
to the synchrotron remnant, or that it is not the 
remnant of SN1181. The ram pressure required 
suggests that the external medium have a number 
density, n ex t, of 

n cxt ~ 0.1 cm~ 3 [u shock /1000km s -1 ] 2 [m/m H ] 

(3) 

where w s hock is difference in velocity between the 
synchrotron bubble and the undisturbed exterior 
medium, and m is the mean mass per particle. 

A possible scenario would be that 3C58 is the 
composite remnant of a binary where both stars 
have undergone supernova explosions, with the 
earlier supernova leaving the slowly expanding 
remnant, and the later, low energy one providing 
the historical supernova event in 1181 A.D. This 
scenario seems rather contrived, since having two 
supernovae within a few thousand years of one an- 
other is quite unlikely 

We will turn to the non-detection of any shell 
emission. If there was an envelope ejected in 



SN1181, as discussed above, it must have been de- 
celerated by at least a factor of 3, which in turn 
suggests that it has already swept up considerable 
external material, and is well into the Sedov phase 
of its evolution. It is then surprising that the ra- 
dio brightness of this exterior shell is considerably 
less then of the faintest known shells. Again, it 
is instructive to point out the difference between 
3C58 and the Crab Nebula. Also for the Crab, 
radio emission from the shell has not yet been de- 
tected to similar limits in surface brightness, but 
in the case of the Crab, the envelope is still in the 
free expansion phase, and is not expected to have 
converted much of its energy into relativistic par- 
ticles, and therefore is expected to be fainter than 
most shell remnants. 

Could a shell for 3C58 be masked because it 
is largely coincident with the synchrotron nebula? 
Shock accelerated radio-synchrotron emitting elec- 
trons typically produce radio emission with a spec- 
tral index, a, ~ —0.6. Our spectral index results, 
however, indicate that there is no significant ra- 
dio emission with a = —0.6 near the edge of the 
nebula, since we observe no steepening of the spec- 
trum there. Also, the total flux density we mea- 
sure at 74 MHz gives an integrated spectral index 
no steeper than that at higher frequencies, con- 
trary to what would be expected if there were a 
steep spectrum component present. 

Both Q!q q7 and show no spatial varia- 

tion over the nebula, ag § 3 showing rms variations 
of only < 0.05. The uniformity of the spectral 
index suggests that a single acceleration mecha- 
nism is responsible for all the radio emitting elec- 
trons. This is thought to be an as yet undetected 
pulsar, although there arc a number of problems 
with 3C58 in this regard. Recent X-ray results by 
Bocchino et al. (2001) do not show the expected 
black-body X-ray emission from the neutron star. 
No pulsed emission has so far been found at any 
wavelength, although in a deep radio image, Frail 
& Moffett (1993) detected a "wisp" like feature, 
similar to ones usually associated with the termi- 
nation shock of a pulsar wind. 

6. A SHELL AROUND 3C58 

In case of the Crab Nebula, there is substan- 
tially less kinetic energy in the remnant than the 
10 51 ergs expected from a type II supernova. This 
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has suggested the existence of a rapidly moving 
shell around the presently visible remnant which 
carries the excess energy and also confines the rel- 
atively high-pressure synchrotron bubble. Much 
observational effort has been expended in search 
of this shell around the Crab, but no shell has 
yet been detected in the radio, ruling out any ra- 
dio emission with a surface brightness much less 
than that of the faintest known young shell rem- 
nant, SN1006 (Frail et al. 1995 and Bietenholz et 
al. 1997). However, in the case of the Crab, cir- 
cumstantial evidence for the existence of expand- 
ing material outside the presently visible nebula 
has been made from Hubble Space Telescope ob- 
servations (Hester et al. 1996) which show the 
Rayleigh- Taylor unstable interface between the 
high-pressure synchrotron bubble and the hydro- 
gen shell. 

In the case of 3C58, there is a similar deficit in 
energy, especially in view of the low expansion ve- 
locity of both the synchrotron bubble and the opti- 
cal filaments, and so a rapidly moving shell might 
also be expected exterior to the presently visible 
nebula. Several searches for a radio shell around 
3C58 have been made, notably by Reynolds & 
Aller (1985), and no sign of radio emission from a 
shell has yet been detected. 

Our 327 MHz results give a sensitive upper 
limit on the radio emission from a putative shell 
of 1 x 10~ 21 W m -2 Hz -1 sr _1 at that frequency. 
Scaling with an assumed shell spectral index of 
—0.6 this is equivalent to 4 x 10~ 22 W m~ 2 Hz -1 
sr _1 at 1.4 GHz, thus very slightly more sensitive 
than the earlier limit of 4.7 x 10~ 22 W m~ 2 Hz -1 
sr _1 obtained by Reynolds & Aller (1985). We are, 
however, more sensitive to steep-spectrum shells 
and also to larger shells due to our larger field- 
of-view at 327 MHz. We note that, even on our 
74 MHz images, where the primary beam radius 
is <~ 5° , there is no sign of shell emission, albeit at 
a lower sensitivity than at 327 MHz. A sensitivity 
to larger shells may be germane if 3C58 is, in fact, 
older than SN1181 since at 10,000 km s _1 , a 3500 
year old shell at 3.2 kpc would have an angular 
diameter of 1.2°. 

7. CONCLUSIONS 

New observations at 74 and 327 MHz and re- 
analysis of older observations dating back to 1974 



indicate that: 

1. 3C58 has a flat spectral index of o^gy — —0.05, 
with little variation over the nebula. 

2. The expansion of the synchrotron nebula is 
much slower than expected for an explosion date in 
1181 A.D. with constant, undecelerated expansion 
since then. A possible explanation is that 3C58 is 
not the remnant of SN1181 but is far older with an 
age of ~ 5000 years, however we are reluctant to 
abandon the generally accepted association until a 
better understanding of the dynamics of plerions 
is available and a better candidate remnant for 
SN1181 is found. An alternate explanation is that 
3C58 has been greatly decelerated. 

3. We are also able to show that, as in the Crab 
Nebula, any non-thermal emission a from a shock 
wave outside the observable plerion has very low 
surface brightness of < 1 x 10~ 22 W m~ 2 Hz -1 
sr _1 . If 3C58 has been greatly decelerated as sug- 
gested in 2. above, then one would expect radio 
emission from such a shock wave. 

We will close by summarizing several appar- 
ent inconsistencies in the observational material 
on 3C58. First, the measured low proper motions, 
both in the radio and in the optical, along with the 
low synchrotron break frequency suggest an age 
much older than 820 years. On the other hand, 
the positional association with SN1181 seems se- 
cure, and no other remnant has been identified for 
this historical supernova. The radio spectral in- 
dex is very uniform, suggesting a common source 
for the radio-emitting electrons, namely a pulsar 
wind, yet the postulated pulsar has escaped detec- 
tion and the expected X-ray emission from the hot 
surface of the neutron star has not been detected. 
Finally, the presently observable remnant contains 
only a small fraction of the kinetic energy expected 
from a type II supernova. Unless SN1181 was an 
anomalously low energy event, a shell consisting 
of the rapidly moving external layers of the star 
would be expected to carry the balance of the ki- 
netic energy. Such a shell has not been detected, 
and is anomalously faint in comparison to other 
known shells. Future study will hopefully allow 
us a more complete and consistent picture of this 
fascinating object. 
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Table 1 

Images used for the Expansion Determination 



Date 




Telescope 


Freq. 


Stotal 3, 


Image rms 


Convolving Beam 


Reference 








(MHz) 


(Jy) 


(mjy/bm) 




°) 




Feb.-Jun. 


1973 


WSRT 


4995 


30.1 


2.44 


8.3 x 8.3 




Wilson & Weiler 1976 


15 Dec. 


1973 


WSRT 


1415 


33.1 


1.35 


22.7 x 20.9 (S 


! 0.0 


Wilson & Weiler 1976 


Jan.-Aug. 


1984 


VLA 


4866 


29.5 


0.43 


8.3 x 8.3 




Reynolds & Aller 1988 






11 


11 


11 


1.71 


17.7 x 17.4 (5 


! -7.5 


11 


Jan.-Dec. 


1984 


VLA 


1446 


33.9 


0.49 


17.7 x 17.4 (Q 


! -7.5 


Reynolds & Aller 1988 


1 Oct. 


1998 


VLA 


327 


33.0 


1.17 


17.7 x 17.4 (5 


! -7.5 


this paper; B- array only 


Mar.-Oct. 


1998 


VLA 


327 


34.0 


0.58 


8.3 x 8.3 




this paper; A+B arrays 



a S to tai is the total cleaned flux density in the image. 
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Table 2 
Expansion of 3C58 



First Image 


Second Image Resolution 81 


Time Interval 


Change in Size 


Rms of Fit 


Weight b 


Year c 


Freq. 


Year c 


Freq. 












(MHz) 




(MHz) (") 


(years) 


(%) 


(mjy/bm) 




1984.39 


1.4 


1998.75 


0.3 17 


14.4 


0.80 


1.4 


1.1 


1984.27 


4.9 


1998.46 


0.3 8 


14.2 


0.45 


0.6 


1.7 


1973.96 


1.4 


1998.75 


0.3 21 


24.8 


0.80 


3.0 


0.5 


1973.96 


1.4 


1984.39 


1.4 21 


10.4 


-0.10 


2.2 


0.9 


1973.96 d 


1.4 


1984.39 d 


1.4 21 


10.4 


0.05 


1.8 


1.3 


1973.29 d 


4.9 


1998.46 d 


0.3 8 


25.2 


0.15 


0.9 


0.8 


1973. 29 d 


4.9 


1984.27 d 


4.9 8 


11.0 


-0.10 


0.7 


1.3 


The results below ; 


ire included 


as a check 










1984.27 


4.9 


1984.39 


1.4 17 


0.1 


-0.20 


1.4 


1.4 


1973.27° 


4.9 


1973.27° 


4.9 8 


0.0 


0.20 


0.9 


0.8 



a The rounded average resolution, see Table 1. 

b These are the weights used for fitting the overall expansion; see text for details. 

°For observations involving multiple arrays, the date given is the average over all the arrays, weighted by 
the length of the individual array observing runs. 

d These images have been high-pass filtered using a circular Gaussian of FWHM 60". 

e This line pertains to the difference between CLEAN and MEM images from the same data 



11 



200 400 




02 02 45 30 15 00 01 45 30 15 00 

RIGHT ASCENSION (B1950) 



Fig. 1. — Image of 3C58 at 74 MHz from combined VLA A & B-array data. The FWHM resolution was 
26" and is indicated at lower left, the peak brightness was 0.55 Jy bm -1 , and the background rms was 
20 mJy bm -1 . The greyscale is in mJy bm _1 . The contours are drawn at —6, 6, 10, 14.1, 20, 28.2, and 40% 
of the peak brightness. 
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Fig. 2.— Image of 3C58 at 327 MHz from VLA A and B-array data. The FWHM resolution was 8.2" and is 
indicated at lower left, the peak brightness was 63 mJy bm _1 , and the background rms was 0.57 mJy bm _1 . 
The greyscale is in mJy bm _1 . The contours are drawn at —2, 2, 2.8, 4, 5.7, 8, 11.3 and 16 mJy bm _1 . 
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Fig. 3. — A low-resolution image of the region surrounding 3C58 at 327 MHz. The image has been convolved 
with a Gaussian of 3.3' FWHM, corrected for the primary beam response, and blanked where the primary 
beam response is < 50%. The peak flux in the image was 15.9 Jy bin -1 . The background rms ranged 
from ~ 20 mjy bin -1 near 3C58 to ~ 45 mJy bm" 1 near the edge of the field of view. The greyscalc is in 
mjy bm^ 1 . 
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Fig. 4. — The radio spectral index of 3C58 between 74 and 327 MHz, afro?, shown in greyscale. The FWHM 
resolution is 21" (shown at lower left). For reference, we repeat the total intensity contours from Figure 1. 
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Fig. 5. — The radio spectral index of 3C58 between 327 MHz and 4.9 GHz, ao'33, shown in greyscale. The 
FWHM resolution is 8.3" (shown at lower left). For reference the 4, 8 and 16% contours from the total 
intensity image at 327 MHz are shown also. 
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Fig. 6. — The expansion of 3C58 over various time intervals. The points plotted represent the expansion 
required by IMDIFF over the pairs of observations from Table 2, excluding the last two check points. 
The solid line indicates the fit expansion speed of 0.020 (±0.08)% yr _1 , and the dashed line indicates the 
expansion of 0.124% yr -1 expected for undecelerated expansion since 1181 A.D. 
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Fig. 7. — Profiles through 3C58. The top panel shows the profile through 3C58, drawn at p. a. —86°, ie. 
through the long axis of the nebula. The solid line represents the flux density at 1.4 GHz in 1984. The 
dashed line represents the flux density at 0.3 GHz in 1998, multiplied by a scaling factor so as to most closely 
match the 1984 profile (the scaling factor was determined by IMDIFF, see text § 4). The dotted line, drawn 
for comparison, is the 1984 data homologously expanded by 1.74%, which is the expected amount if 3C58 is 
the undecelerated remnant of SN1181. The middle and the lower panel are magnified regions of the profile in 
the top panel, showing details of the eastern and western limbs of 3C58 respectively. The FWHM resolution 
was 17.7" x 17.4" @ -7.5°. 18 



